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ABSTRACT 

This paper presents a simplified model of the space shuttle main engine (SSME) 
dynamics valid within the range of operation of the engine. This model is obtained 
by linking the linearized point models obtained at 25 different operating points of 
the SSME. The simplified model was developed for use with a model-based diagnostic 
scheme for failure detection and diagnostics studies, as well as control design 
purposes . 


INTRODUCTION 

This paper presents a simplified mathematical model describing the dynamics of 
the space shuttle main engine (SSME). A complete nonlinear dynamic simulation of the 
SSME was developed by Rockwell International Corporation (ref. 1). Due to its size 
and complexity, however, this nonlinear simulation is difficult to use for various 
design and diagnostic purposes. Therefore, a simplified dynamic model of the SSME, 
valid within the range of engine operation, is required for model based diagnostic 
and control system design. 

The description of the SSME and its control system is presented in a work by 
Duyar, Guo, and Merrill (ref. 2). Also in a previous study, Duyar, Eldem, Merrill, 
and Guo (ref. 3), developed an identification algorithm and used it to obtain 
linearized models of the SSME at five different operating points. There the five 
different operating points were obtained over a range of main combustion chamber 
pressures at a constant value of mixture ratio. This study is a continuation which 
increases the number of operating point models to 25 and links perturbation models 
identified at these operating points to obtain a simplified model of the SSME. Here, 
the operating points not only include a range of main combustion chamber pressures 
but also a range of mixture ratios. The coefficient matrices of the linear models 
were regressed with the operating point parameters (chamber pressure and mixture 
ratio) to yield a simplified nonlinear model valid within the range considered. The 
accuracy of the simplified model was tested by designing an input signal which covers 
the whole range of operating points considered and then comparing the responses of 
the linked model and the nonlinear simulation. This approach has been successfully 
applied to control and diagnostic design in other propulsion systems (ref. 4). 

Possible applications of this model include a model-based diagnostics scheme for 
failure detection purposes which compares the parameters of the simplified model with 
the identified parameters of the actual system on a continuous basis. The 
differences between the modeled and identified parameters may be analyzed to detect 
the presence, to isolate the source, and to determine the extent of SSME component 


failures. Other applications for this model would include predictions of remaining 
component life, open loop engine dynamics studies, and closed loop control analysis 
with a user generated control law. 

The remainder of this paper includes a brief description of the data acquired 
from the nonlinear dynamic simulation. This is followed by a description of the 
structure of the identified point models and the linked model. Finally, simulated 
results obtained from the linked model are compared with the results obtained from 
the nonlinear simulation for the same test input. 


DATA ACQUISITION 

In order to obtain input-output data of the SSME, a nonlinear dynamic simulation 
program was used. The main combustion chamber mixture ratio, MR n , and pressure, P cn , 
and the engine face pressure, P , are the parameters which define a nominal operating 
condition. However, it was observed that the engine face pressure does not signifi- 
cantly affect the selected state of the output variables of the system. Hence, only 
the mixture ratio and the chamber pressure were selected as the defining variables. 

To obtain the data at a particular operating point, the nonlinear simulation 
with closed-loop control was allowed to converge to the desired operating point. 

Then, the control was removed and the actuator outputs, oxidizer preburner oxidizer 
valve (OPOV) rotary motion and fuel preburner oxidizer valve (FPOV) rotary motion 
were disturbed simultaneously using two uncorrelated full length, pseudo-random 
binary sequences (PRBS), The maximum magnitude of each sequence was 2 percent of the 
steady-state value of the input. The clock-time interval chosen for each sequence 
was 0.04 sec. The length of the PRBS used was 127 intervals for each input which 
corresponds to a signal duration of 127 x 0.04 = 5.08 sec. Thus, the chosen clock- 
time and signal duration give a sequence with a bandwidth from 0.2 to 12.5 Hz. 


POINT MODELS 

Once the input/output data (of length 127) was obtained at each of the operating 
points, the developed identification algorithm was used to obtain the point models. 
First, a two-input/two-output model was considered. Here the outputs were the devi- 
ations of the chamber inlet pressure, P c , and the mixture ratio, MR and the inputs 

were the deviations of the valve actuator output rotary motion B and B 

r i r opov ~ FPO v 

From the operating point data the identification algorithm yielded a fourth order 
system in a-canonical form with observability indices (2,2). This result encouraged 
further investigation to obtain a model with measurable state variables. To this end 
the linear dependences between several different outputs and P and MR were 
compared and finally high pressure fuel turbine, HPFT, and high pressure oxidizer 
turbine, HPOT, speeds, S and S , were chosen as third and fourth state 

HPFT HPOT 

variables for the dynamic model of the SSME. Then, the identification algorithm was 
employed to determine the parameters of this conceptual model. This model can be 
expressed mathematically as: 
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( 1 ) 


Sx^k +1) = A i <5x t (k) + 5u l (k) 


i = 1,2,. ..,25 


where A and B i are 4x4 and 4x2 system matrices at operating point i and 
<5x < . ) is the model estimate of the state, which is also the model estimate of the 
output, since C = I. The perturbed input is defined as 


5u iV° 


A ± j 5u 1;) (. ) , 5u i;) (. ) l o 

(2 — ( • ) r ^ 0 


(2) 


with 5u li (.) and 5u ±2 ( . ) denoting the perturbed input values /? opov and /? FP0V at 
operating point i, respectively. Perturbed values of the inputs are 

<5u i;j (. ) = (U tj (.) - U^j/u^ ( 


where u ± ^ denotes the steady state value of input j at operating point i. In a 
previous study (ref. 3) it was observed from the open loop response of the SSME that 
the gains for positive and negative perturbation signals are significantly different. 
In order to compensate for this phenomenon, piecewise linear system gains, A^, with 
different values for positive and negative perturbations were added to the linearized 
model. These piecewise linear gains defined by equation (2) were determined by using 
the sample averages of outputs and inputs at each operating point as explained in 
reference 3. 

The perturbed values of the measured outputs of the nonlinear simulation, 

5y ( . ) r at each operating point were defined as 

5y n (.) = (y n - Ynl/Yn (4) 


where 

y deviation of the chamber inlet pressure at operating point i 

y steady state value of the chamber inlet pressure at operating point i 

Similar definitions are valid for 6y ( . ) where the index j = 2,3,4 denote 
MR, S and . respectively. 

HP FT HPOT c * 

The parameter matrices A i and B 1 were also employed to obtain a state filter 
estimate tfz^.) of 5y i (.) using the following equation: 

6z i (k +1) =A 1 5y i (k)+B i 5u i (k) i * 1, . . . , 25 ^ 


THE LINKED MODEL 

The coefficients of the 25 operating point models were regressed with the steady 
state values of P cf} and MR n . The variables P cn and MR n define an operating 
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point. A second order (with respect to P cn and MR n ) polynomial expression for 
each coefficient was obtained. These 25 operating points correspond to five MR n 
levels and five power levels at each MR^ level as follows: 


1-5 

correspond to 

MR = 6.026 

N 


and 

P CN ” 

70, 

O'* 

o 

CO 

i, 100, and 105 percent 

6-10 

correspond to 

MR n = 6.026 

X 

1.04 

and 

P CN 

= 70, 

80, 

90, 100, 

and 

105 

percent 

11-15 

correspond to 

MR =6.026 

N 

X 

0.96 

and 

P CN 

= 70, 

80, 

90, 100, 

and 

105 

percent 

16-20 

correspond to 

MR =6.026 

N 

X 

1.08 

and 

P CN 

= 70, 

80, 

90, 100, 

and 

105 

percent 

21-25 

correspond to 

MR n = 6.026 

X 

0.92 

and 

P CN 

= 70, 

80, 

90, 100, 

and 

105 

percent 


In regressing the coefficients however, the steady state values of the inputs 
and outputs at each operating point should be taken into account. This can be under- 
stood by substituting the perturbation expression of equations (2) and (3) into equa- 
tion (1). Note that (Sx^k) is the deviation of the model estimate of the state 
vector at the i th operating point, normalized with respect to the steady state 
values at that operating point defined by component as 




Xjfk) - x i;j 


X 


ij 


( 6 ) 


where the subscript j denotes the j th component of the state vector and x 1 . is 
the steady state value of the j th component of the state at the i th operating 
condition. Thus equation (1) can be rewritten as 


E 


* x (k + 1) 



A i E i 


x(k) 




(7) 


where x ± and u ± are the steady state values of the outputs and inputs of the i th 


operating point. Here, x ± denotes the estimate of y 1 . 
as 

E i and T i are defined 

E i = diag {x i;) } 

(8) 

r i = dia< 3 {^i]} 

(9) 


Further rearrangement of equation (7) gives 
x ( k + 1) = E 1 A 1 E’ 1 x(k) + E^r" 1 

Thus, the coefficients to be regressed for the linked model are: 

E i A 1 E i 1 , EJjJ’/, X ijf u lf and x 1 . 

The regression yields the linked model as 


u(k) - Ul (k) 


* + (I - 


( 10 ) 
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x(k + 1) = A[x(k)]x(k) + B[x(k)]{u(k) - u[x(k) ]}* + {i - A[x(k) ]}x[x(k) ] 
for the model, and similarly 

z(k + 1) = A[y(k) ]y(k) + B[y(k) ] {u(k) - u[y (k) ] } * + {i - A[y(k> ]}y[x(k) ] 


for the linked filter. In the above equations 


u. ( k) - u,( . ) 


Xjfofu^k) - u^.) , 

[2 - \)(0][uj(k) - u^.) 


u d (k) l u d ( . ) 

, u d (k) $ u d ( . ) 


(13) 


where (.) denotes a function of x(k) for the model and y(k) for the filter. 
Furthermore, 


A[ . ] * A c + A„P 


P‘ C + ^ + ^p2^C + ^M2 ^ + ^PM^C 


MR 


(14) 


where P and MR are the current (k th value) of the main combustion chamber inlet 

c 

pressure and the mixture ratio. For the model the estimated values at the previous 
stage were used, i.e., A[ . ] = A[x(k)]. For the filter the actual outputs from the 
nonlinear simulation were used, i.e., A[ . ] = A[y(k)]. Note that A c denotes a 

constant term. Similar expressions can be written for B, u, and x which yield: 


B c' B p' V B P2' B M2' B PM f Sc' *jp' ^JP2' ^ jM2 ' \m' 


for j a 1,2 (for two inputs) 


The above coefficients were obtained by regressing the coefficients of the linear 
models with P cn and MR n values of each operating point and are presented in 
table 1. 

The regression was performed using the Statistical Package for Social Sciences 
(SPSS) (ref. 5). The chosen second degree polynomial relationship is shown in 
equation (14). The objective here was to include only those independent variables 
that were useful in predicting the dependent variable. The partial F test 
criterion was used to determine the contribution of an independent variable. The F 
statistic was calculated by testing the null and alternative hypothesis as follows: 

H q : There is no relationship between the dependent variable and the independent 

variable under test (the coefficient under test is 0) . 

H : There is a definite relationship between the dependent and the independent 

variable (the coefficient under test is not 0). 

Using a 95 percent level of significance, the null hypothesis is rejected if 
F > F crltical * The stepwise regression procedure in SPSS is used to perform such a 
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search. Note that the equations which define the linked model and the filter have 
linear structure but the coefficients are quadratic functions of MR and P . 

N CN 


MODEL VERIFICATION 

The linked model was tested for accuracy by comparing its response to a test 
input signal with the response obtained from the nonlinear simulation. 


The input test signal was designed such that it covers the whole range of 
operating points considered. For this purpose, input test signals for P cr and MR r 
requests were designed first. These signals consist of steps and ramps ranging from 
70 to 100 percent power levels and mixture ratio levels of 5.54 to 6.51. Figure 1(a) 
shows the 
and 


MR r and P cr request signals and figure 1(b) shows the corresponding 

B actuator commands generated by the controller. (Note that the same 

ranges for P cr and MR r were used in identifying the point models. ) Then the 
OPOV and FPOV values corresponding to this closed-loop operation shown in 
figure 1(b) are used as inputs to the linked model (for the filter, P c and MR 
values are also used) . For a visual comparison, the outputs of the nonlinear 
simulation, y(k) , the model x(k) and the filter, z(k), are shown in figures 2 and 3. 
Standard error of estimates (SEE) for both the model and the filter were also 
calculated according to the formula 


SEE 


E[y(k) - x ( k ) ] 2 


|l/2 


£[y(k) j 2 


(15) 


Table II gives the SEE's of the predictions of both the linked model and the 
filter. It was observed that both the model and the filter have a good fit to the 
data from the nonlinear simulation. 


CONCLUSION 

An off-line algorithm presented in a previous study (ref. 3) and the data 
generated from the nonlinear performance simulation (ref. 1) were used to obtain 
linearized models of the SSME at 25 different operating points. The models have 
measurable state variables which make them convenient for control design and failure 
detection and diagnostic studies, since the need for state estimation is eliminated. 

These point models were linked to obtain a simplified model of the SSME, valid 
within its full range of operation. The coefficients of the linear models were 
regressed with the parameters determining the nominal operating conditions: the 

chamber mixture ratio, MR n , and pressure, P cN . The engine face pressure, P , was also 
considered but it was observed not to have a significant effect on the selected state 
and the output variables of the system. The regression yields a nonlinear model 
valid within the range of mixture ratio and power levels considered. This simplified 
model was tested by designing an input signal which covers the whole range of 
operating points considered and then comparing the responses of the linked model and 
the nonlinear simulation. The comparison indicated good agreement between these 
responses. It is concluded that the simplified model is sufficiently accurate to be 
used in model-based diagnostics. 
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TABLE 1. - THE PARAMETERS OF THE LINKED MODEL 


Coeffi- 

cients 

A 


B 


X 


X 


0 

A 

1 . 3072 

0.2404 

7.8608 

0.5497 

B 

0.4453 

0.0074 

x r 



0.0000 




0.1604 

1 . 2022 

15.2350 

0,3455 


-21.1840 

0. 1509 


0.9505 

X c 

0.0000 

u c 

0.8922 


0.9131 

0.1321 

0.5400 

0.2032 


0. 1853 

1. 1043 


0.8790 


0.8644 


1.0504 


0.0701 

0.0650 

0.3075 

-2.5536 


-1.6421 

0.1637 




0.1966 



a m 

0.0000 

-0,6780 

-12.2966 

0.0000 

B m 

0.0000 

0.0000 




0.0000 

o 



0.0000 

0.0000 

-27.9998 

1.3465 


38.2470 

0 . 0000 


0.0000 

*M 

1.0000 

U M 

0.0000 


0.0000 

-0.1238 

0.0000 

0.0000 


O.OCTOO 

-1.9433 


0.0000 


0.6326 


0.0000 


0.0000 

0.0000 

0.0000 

5.7173 


0.0000 

0.0000 




0.0000 



A 

-2.4623 

0.3456 

-3.7077 

0.0000 

Bp 

0.0000 

0.0000 

K 


Xp 

1.0000 

Up 


P 

0.0000 

-2.4590 

-5.1573 

0.0000 


8. 3042 

0.0000 


0.0000 


0.0000 


-0.5759 


-2.0943 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 


0.0000 


0.5529 


0.0000 


0.0000 

0.0000 

0.0000 

0.0000 


6.0243 

-0.2370 




0.6092 




1.5139 

0.0000 

0.0000 

0.0000 

B P2 

0.0000 

-0. 1625 

Ki 



0.0000 



P2 

0.0000 

1.3207 

0.0000 

-0.8525 


0.0000 

0.0000 


0.0000 


0.0000 

0.6099 


1 . 1225 

0.0000 

0.0000 

0.0000 


0.0000 

-2.3360 


0.0000 


0.2657 


0.6450 


0.0000 

0.0000 

0.0000 

0.0000 


-3.8452 

0.0000 




0.0000 




0.0000 

0.4078 

4.8210 

0.0000 

B M 2 

0.0000 

0.0000 



x 

0.0000 



M2 

0.0000 

0.0000 

9.2634 

0.0000 


-14.2340 

0.0000 


0.0000 

A H2 

0.0000 

M2 

0.0798 


0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 


0.0000 


0.3665 


0.0000 


-0. 1066 

-0.0437 

0.0000 

-2.7117 


0.0000 

0.0000 




0.0000 



A 

-0.3621 

-0.2947 

3.9062 

0.0000 

B pm 

-0.2277 

0.1386 



7 

0.0000 



PM 

0.0000 

0.0000 

6.8391 

0.0000 


-10.0470 

-0.2312 


0.0000 


0.0000 


0.0000 


0.0000 

0.0000 

0.0000 

0.0000 


-0.1127 

3.4060 


0.0000 


-0.4159 


-0.6671 


0.0000 

0.0000 

0 . 0000 

0.0000 


0.0000 

0.0000 




0.1913 

— 
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Normalized i 


TABLE 2. - STANDARD ERROR OF ESTIMATES OF 
PREDICTIONS OF THE FILTER AND THE MODEL 
FOR THE TEST SIGNAL COVERING ALL 
OPERATING POINTS 


Test signal 

Output 

P 

c 

MR 

S 

HPFT 

S 

HPOT 

Filter 

Model 

0 ♦ 0039 
.0128 

0.0114 

.0244 

0.0076 

.0151 

0.0028 

.0075 




Figure 1 .—Reference input signals and their corresponding valve inputs. 
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Normalized output 


1.10 


Simulation 

Filter 

Model 
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